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Abstract 

Ground state properties of the superheavy elements (SHE) with 108<Z<128 
and 150<iV<192 are investigated using both the Skyrme-Hartree-Fock 
method with a density-independent contact pairing interaction and the 
macroscopic-microscopic approach with an average Woods-Saxon potential 
and a monopole pairing interaction. Detailed analysis of binding energies, 
separation energies, shell effects, single proton and neutron states, equilib- 
rium deformations, Q a -values, and other observables is given. 

PACS number(s): 21.10.-k, 21.30.Fe, 21.60.Jz, 23.60.+e, 27.90.+b 
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I. INTRODUCTION 



Recent years have witnessed great strides in the production of the heaviest nuclei. No- 
tably, three new elements, Z=110, 111, and 112, were synthesized at GSI [l]]^, Berkeley 0, 
and Dubna 0] . The systematics of reaction data show that cross-sections are monotonically 
decreasing, reaching a level of lpb for production of element 112. The half-lives of the 
known isotopes of elements with Z>105 are predominantly determined by the a-decay, and 
decrease from 0.9s for 263 1 06 to 0.2ms for 277 112 165 . These isotopes decay predominantly 
by groups of a particles (a-chains). The corresponding a-spectra carry important nuclear 
structure information. 

All the heaviest elements found recently are believed to be well deformed. Indeed, the 
measured a-decay energies, along with complementary syntheses of new neutron-rich iso- 
topes of elements Z=106 and Z=108 by the Dubna-Livermore collaboration, have furnished 
confirmation of the special stability of the deformed shell at iV=162 PHT5]]. Still heavier and 
more neutron-rich elements are expected to be spherical and even more strongly stabilized 
by shell effects; they form the long predicted region of the superheavy elements (SHE). A 
long-standing goal in themselves, SHE are also a critical test of nuclear models, of relativistic 
effects, and of quantum chemistry. 

Experimentally, the excursion towards the center of the SHE will not be easy. Indeed, 
for the most neutron-rich isotopes synthesized, production cross sections are of the order of 
pb which is on the limit of today's experimental sensitivity. Two ways of synthesizing new 
elements have been employed. Experiments at GSI utilized cold fusion reactions 208 Pb(HI, 
xn) which benefit from the enhanced stability of the target nuclei. On the other hand, the 
Dubna-Livermore experiments utilized more asymmetric hot fusion reactions 244 Pu(HI,xn). 
As far as future measurements are concerned, there have been intensive discussions on how 
to reach still-heavier new nuclei P|16|-|i9|. Though the experiments will be difficult, and 
will probably require advances in separation and detection technology, one advantage is 
that many of the nuclei sought are predicted to have longer lifetimes than those recently 
synthesized. Clearly, there is now a wealth of possibilities opening up involving cold and 
hot fusion, light and heavy projectiles, and stable and radioactive beams. 

Among various theoretical approaches to the heaviest elements, particularly successful 
have been calculations based on the macroscopic-microscopic method (Nilsson-Strutinsky 
approach). In particular, detailed studies have been carried out with the folded- Yukawa (FY) 
deformed single-particle potential [p0| -27l and with the Woods-Saxon (WS) deformed single- 
particle potential [pj- pl]JI^ , ^5|j28| , [29| . In most cases, the FY and WS results are consistent; 
they both reproduce well-known alpha-decay half-lives of heavy elements and emphasize the 
stabilizing role of deformation for nuclei with N~162. Since the liquid-drop barriers in SHE 
disappear, the shell effects dramatically influence the spontaneous fission and alpha-decay 
half-lives; the delay in the spontaneous fission half-lives, T s f, due to the shell effects is about 
15 orders of magnitude for Z>106 |]I6|. In this context, it is important to recall that in both 



FY and WS models the center of the shell stability is predicted around the alleged spherical 
doubly- magic nucleus 298 114 184 . 

The aim of the present study is to understand the shell structure properties of SHE based 
on the self-consistent theory. Although the physics of the heaviest and superheavy nuclei has 
been addressed by numerous calculations (for a review, see Refs. @,0), a large-scale self- 



2 



consistent description of SHE involving non-spherical degrees of freedom is still lacking. It is 
our intention to fill this gap. The macroscopic-microscopic approach, which has mostly been 
used up to now, is only an approximation to the Hartree-Fock (HF) theory and suffers from 
the lack of self-consistent coupling between the macroscopic and microscopic contributions 
to the total energy. Moreover, it is not based on a microscopic many-body Hamiltonian and 
its success (or failure) in reproducing experimental data cannot be easily traced back to an 
effective interaction (or to the corresponding density functional). In this work we analyze 
the properties of the SHE by means of the self-consistent Skyrme-Hartree-Fock approach 
with a zero-range pairing interaction. For comparison, a set of calculations has also been 
carried out with the macroscopic-microscopic WS approach. 

This work has several objectives. Firstly, we re-examine the old problem of the shell 
stability of SHE. Namely, we study single-particle spectra resulting from the (unusually 
strong) self-consistent interplay between the nuclear and Coulomb parts of the Hamiltonian, 
discuss the major spherical gaps, and calculate ground-state shapes. Secondly, we analyze 
the impact of the self-consistent treatment on lifetimes of SHE and, in selected cases, make 
estimates of fission barriers. Last but not least, we systematize predicted binding energies, 
separation energies, and a-decay half-lives and discuss their dependence on the choice of the 
Skyrme parametrization. 

The paper is organized as follows. The models employed are described in Sec. ||. Section 
|H contains the analysis of the single-particle structure of SHE. Ground-state quadrupole and 
hexadecapole deformations and shell corrections are discussed in Sees. [TV] and [V], respectively. 
Binding-energy differences (i.e., separation energies and Q a values) are analyzed in Sec. [\T|. 
Section |V11| deals with the particle-decay lifetimes while Sec. |V111| contains selected results 



of fission barrier calculations. Finally, conclusions are contained in Sec. IX 



II. THE MODELS 

The calculations presented in this paper are based on the mean-field theory. The self- 
consistent results have been obtained using the Skyrme-Hartree-Fock+BCS method and the 
macroscopic-microscopic results with the deformed WS model. These models are described 
in Sees. |II A| and |IIB| , respectively. 



A. Hartree-Fock calculations 



The Hartree-Fock equations have been solved according to the method of Ref. |JT]. The 



main advantage of this method is that instead of expanding single-particle wave functions on 
an oscillator basis, one discretizes them on a three-dimensional mesh, with equally spaced 
points along three cartesian axes. This method has been shown to have an accuracy nearly 
independent of the nuclear deformation, and it does not require an additional basis opti- 
mization. The high efficiency of the method can be attributed to the iterative procedure 
used to solve the HF equations, the imaginary time step method. One of the features of 
this method is that only those single-particle wave functions which contribute to the mean 
field are calculated. In practice, in the following we have considered all the wave functions 
which have a BCS occupation larger than 10~ 5 , i.e., approximately 160 single-particle states 
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for the protons and 280 for the neutrons. The actual mesh consisted of 15 3 points with a 
mesh step of 1.0 fm. Under such conditions, the accuracy of calculations can be estimated to 
be equivalent to that obtained with an oscillator basis expansion in optimized 18 deformed 
shells. 

In the calculations we imposed the existence of three symmetry planes. The resulting 
shapes are reflection-symmetric, but they can be triaxial. These assumed self-consistent 
symmetries are expected to be well justified for deformations slightly beyond the second 
minimum. For larger elongations, however, one expects the presence of odd-multipolarity 
moments, disregarded in this work. 

In the particle-hole (p-h) channel, the Skyrme effective interactions have been used. The 
two parametrizations employed here, SkP and SLy7 [|34l,|35| , satisfactorily describe 



the properties of exotic nuclei far from the beta stability line. Moreover, these interactions 
have been adjusted to reproduce long isotopic sequences; hence one can expect them to 
have good isospin properties. SkP has been especially designed to act both in the p-h and 
particle-particle (p-p) channels. It has an effective mass equal to the free nucleon mass 
(m*/m=l). Therefore, it reproduces properly the single-particle level density around the 
Fermi surface. When fitting SLy7, special attention was paid to its ability to describe the 
properties of symmetric and asymmetric infinite nuclear matter. For this force the effective 
mass is rather low, m*/m=0.67. Surface properties of SLy7 are close to those of the SkM* 
parametrization. Consequently, SLy7 is expected to describe correctly heavy nuclei at large 
deformations, in particular fission barriers. During the adjustment of the parametrization 
SLy7, both the tensor contribution to the spin-orbit term and the two-body center-of-mass 
correction have been included self-consistently on the mean-field level. In our study, however, 
these terms are treated perturbatively, i.e., the corresponding corrections to the mean fields 
are not taken into account and only added to the total energy after the self-consistent 
solutions are found. Such a procedure may slightly affect the single-particle levels obtained 
in the HF+SLy7 variant discussed in Sec. pTT| . 

The nuclear matter characteristics of SkP and SLy7 are displayed in Table |. These forces 
differ mainly by different values of the isoscalar effective mass and enhancement factor (i.e., 
isovector effective mass). The remaining nuclear matter properties are very similar. In par- 
ticular, the corresponding semiclassical surface energies |$6f are almost identical and equal 



to that for the SkM* interaction ]37| a sur f=21.8 MeV, where it has been adjusted to give 



correct fission barrier properties. More refined methods to calculate a sur f give slightly 
different numbers which, however, are again very similar for the three forces SkM*, SkP, 
and SLy7. Nevertheless, the fission barriers obtained for SkP and SLy7 are significantly 



different, Sec. IVTTJ. This suggests that the nuclear matter indicators are too crude a handle 



on the results of actual HF calculations at large deformations where a detailed shell struc- 
ture is more important than the overall semiclassical properties. By using the two Skyrme 
parametrizations, SkP and SLy7, which differ in their effective mass properties, we aim at 
investigating the uncertainity in predictions related to the p-h interaction channel. 
In the p-p channel, we use a density-independent contact interaction, 

V T 5 (r,r') =V 0>T 6(r-r') (r = n,p), (1) 

for both Skryme parametrizations employed. In order to minimize unphysical fluctuations of 
the pairing field in regions of low level density, an approximate particle number projection 
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was implemented by means of the Lipkin-Nogami (LN) method |39|.|40| . 

In order to adjust the pairing strengths V 0>T , we have estimated an average A-dependence 
of the experimental pairing gaps for even-even nuclei around 254 Fm. As a result, the following 
expressions for the proton and neutron average pairing gaps have been obtained: 

A r = St/A 1 / 2 , (2) 

where 5 n = 10.304 Me V and <5 P =13.072 MeV. The pairing strengths of the interaction (|]) have 
then been adjusted to reproduce the average proton and neutron pairing gaps in 254 Fm. 

Table [II] displays the resulting pairing strengths, together with the predicted and ex- 
perimental pairing gaps in 254 Fm. The experimental pairing gaps, A exp , were extracted 
from odd-even mass differences taken from Ref. JO]. The quantity A LN + A2 represents the 



"effective" pairing gap in the LN approximation (see, e.g., Ref. [f|2[]). The rather different 
values of V obtained in SkP and SLy7 reflect the differences in their respective effective 
masses. Indeed, the larger values of Vq in SLy7 compensate for the reduced single-particle 
level density around the Fermi level due to the low value of m* ||43|| . As shown in Table |TI], 
the optimized values of Vq )T reproduce rather well the average pairing gaps in even-even 
nuclei in the Th-Fm region. 

It is well known |3^Q] that traditional models of pairing correlations (e.g., BCS) become 
inappropriate when approaching particle drip lines. The main drawback is the scattering 
of nucleons from bound shell-model orbitals to unbound states, giving rise to an unphysical 
component in the nucleonic density with the wrong asymptotic behavior. This problem 
becomes particularly severe for neutron-rich nuclei where the particle continuum lies very 
low in energy. On the proton-rich side, which is approached in the present calculations, the 



effect of the continuum is weaker [H5| because of the Coulomb barrier which tends to localize 



the proton density in the nuclear interior. 



B. Woods-Saxon calculations 



The macroscopic-microscopic method is an approximation to the HF approach [46,4? . 
Its main assumption is that the total energy of a nucleus can be decomposed in two parts, 



E = E n 



(3) 



where -E ma cro is the macroscopic energy and -Emicro is the microscopic energy (shell correction) 
calculated from a non-self-consistent average deformed potential. 

In our study, -Emacro was calculated using the Yukawa-plus-exponential mass formula 
of Ref. [24| with parameters of Ref. | 13[| . That is, similar to Ref. we adjusted the 
volume-asymmetry parameter, k v , and the charge-asymmetry parameter, c a , to obtain good 
local agreement with known masses in the heaviest element region. The optimal values, 
k v =1. 92552 and c a =0. 14505 MeV, turned out to be close to the original values of Ref. p4| . 
i.e., k„=1.911 and c a =0.145 MeV. 

The shell correction -E m i cro was computed using the axially-deformed single-particle 
Woods-Saxon (WS) Hamiltonian ||8],f49| with the parameters of Ref. [ 50(1 - In the WS calcu- 



lations the nuclear surface is defined by means of standard fix axial deformation parameters: 
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R(Q) = c(p)R 



A=2,4,6 



(4) 



with c(/3) being determined from the volume-conservation condition and R = r A 1//3 . The 



shell correction was calculated according to the prescription given in Ref. [J47|, assuming the 



smoothing parameter of 49.2/A 1 / 3 MeV and a sixth-order curvature correction. All single- 
particle states lying below a cut-off energy of 131. 2/A 1 / 3 MeV above the Fermi level were 
included. 

The pairing energy was calculated in the LN method using a state-independent monopole- 
pairing interaction. The pairing strengths and the average pairing energy were taken accord- 
ing to Ref. . The inclusion of the average LN pairing energy led to a renormalization of 
the mass-independent constant in the macroscopic energy formula [2"0|,2"3|, Co=l. 59741 MeV. 
Finally, the zero-point energy correction of Ref. [29] was added to the total binding energy. 



III. SINGLE-PARTICLE STRUCTURE OF SHE 

The shell structure of SHE is governed by several factors, the most important one being 
the large mass and large electric charge of a superheavy nucleus. With increasing A, the 
single-particle spectrum becomes more compressed, and the binding energy increases due 
to the increased range (radius) of the average potential. Moreover, the orbitals with £—J\f, 
j=£—l/2, such as lh 9 / 2 , IJ13/2, etc., are considerably shifted down in energy due to 

a weakened spin-orbit splitting at large values of A (see Refs. [^l],^2j and Fig. 2-30 of Ref. 

Indeed, according to the estimates of Refs. the one-body coupling constant 

of the spin-orbit term decreases faster than expected from the usual radial A~ 1 / 3 -scaling 
law flSql . For protons, an additional change in the shell structure comes from the Coulomb 
potential which decreases the binding energy of the single-proton orbitals with small values 
of (r 2 ). Consequently, the Coulomb potential gives rise to a lowering of the unique-parity 
shell (£=Af) with respect to the normal-parity orbitals with J\f'=J\f—l. 

Some of the systematic changes in the single-particle shell structure are illustrated in 
Figs. ^ and |2| which display the spherical WS single-proton energies in 208 Pb and 310 126i84, 
respectively, as functions of the Coulomb potential strength xcoui- The latter parameter, 
varied in the range 0<xcoui<l ; renormalizes the electric charge Z—l that defines the strength 
of the Coulomb potential. Namely, for £coui=0 the Coulomb term vanishes, and for a;coui=l 
the Coulomb strength takes its normal value. In order not to expand the energy scale too 
much, the single-particle levels are shifted in Figs. p] and |] by 

Ae = ^ C oui(r = 0) [1 - x Co ui] , (5) 

where Vcoui( r = 0) is the Coulomb energy at the center of the nucleus. 

The influence of the Coulomb term is apparent in Fig. |l|. With increasing xcoui, high-j 
shells (most notably: unique parity shells) are clearly shifted down with respect to low-j 
orbitals. This shift, however, does not lead to any significant change in the spherical shell 
structure. In particular, the Coulomb term does not change the position of magic gaps 50, 
82, and 126. The situation is different in the superheavy nucleus 310 126 (Fig. 0). Here, the 
Coulomb energy acting together with the energy shifts due to (i) the increased radius and 
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(ii) the change in the spin-orbit splitting discussed above, induces significant changes in the 
shell structure. In particular, the lowering of the I/19/2 and lin/2 orbitals gives rise to the 
(almost complete) closing of the spherical gaps at Z=82 and 126, and the appearance of the 
subshell closure at Z=114 (cf. also Fig. 2 of Ref. |52| and discussion in Ref. [j56| ). 

The position of the magic gaps in the superheavy region has been discussed extensively 
in the past, both in non-self-consistent [|TT1 , p6| , ^^ , |5T| , |52| , ^5H 7^1 and self-consistent [f75|-|96H 
models. Early calculations based on finite average potentials |^T|J^] and the modified os- 



cillator potential [0,0 found Z=11A (but not Z—126) and iV=184 to be the magic num- 
bers in the superheavy region. (Similar conclusions were reached in the following works 
56| , |65|j67| , [74]j97| . ) The magic character of the Z—114 gap was questioned by several authors 



IBTfl who pointed out that its size is sensitive to the details of extrapolation of model param- 
eters; in particular it is influenced by the magnitude of the 2/7/2 — 2/ 5 / 2 spin-orbit splitting 
|66| , |98|| . The size of the 2=114 and iV=184 gaps in the nucleus 298 114, predicted in various 



models, is displayed in Table [IV]. In almost all cases, the size of the Z=114 gap reflects the 
2/7/2-2/5/2 proton spin-orbit splitting. In some calculations, e.g., HF+SLy7, HF+SkP, and 
Briickner-Hartree-Fock (BHF), the 1*13/2 proton orbital appears above the 2/7/2 shell. (In 
several other models, including WS, HF+SIII these two shells are practically degenerate, 
see Fig. 0; we also note the similarity between the HF spectra of Refs. [BILK8[[52] and our 
HF+SkP spectrum.) 

Generally, calculated shell effects in the SHE (hence their lifetimes) depend very strongly 
on the symmetry dependence [(N—Z)/ A - term] of the spin-orbit potential [^,^]. In this 
context, it is interesting to compare in Table [IV| the Skyrme-HF results with those obtained 
in the RMF theory p4| , p5| , where the the isospin dependence of the spin-orbit term is 



different. 

Figures ^|-^] display the spherical single-particle levels in 294 114 and 310 126 obtained in 
the HF+SkP, HF+Sly7, and WS calculations. For the neutrons, all variants give a rather 
consistent scenario: the first significant spherical shell gap above iV=126 is 184, and the 
next large spherical closure is predicted at iV=228 (see Fig. ||). It is to be noted that the 
neutron shell structure, including the absolute shell-gap sizes, is extremely similar in the 
HF+SkP and WS variants, while the spectrum obtained in HF+SLy7 shows a characteristic 
stretching. This can be attributed to a fairly small effective mass in SLy7; a smaller effective 
mass gives rise to a reduced single-particle level density, hence to an overestimation of shell 
effects |33|j90|j93[1 . For this reason, the SkP spectrum (effective mass equal to one) is close to 
the WS spectrum. 

For the protons (Figs. |] and |^) the results show strong model dependence. Namely, 
the HF+SkP variant yields the magic proton gap at Z=126; in the HF+SLy7 variant the 
sizes of the Z—114 and Z=12Q gaps are comparable (although the level density around the 
Z=126 gap is much lower), and there is a strong preference for the Z—114 gap in the WS 
calculations. Again, the SkP interaction predicts smaller gaps than SLy7. 

In light of a fairly good agreement for neutrons, the visible discrepancy between the 
WS results and the self-consistent results for protons might look surprising. This difference, 
however, has a simple explanation in terms of surface properties. Indeed, it has been argued 
for some time that a serious source of uncertainty, when extrapolating the parameters of 
average potentials such as WS to the region of SHE, is the surface diffuseness or surface 
thickness [B2l,|53|,|56l,|72l,|8^,|S3,S9fl. Since self-consistent models indicate that the surface thick- 
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ness parameter for protons is fairly close to that for neutrons (cf. discussion in Refs. |93| . |10C 



based on the semiclassical extended Thomas- Fermi (ETF) approach), it has been argued 
|7^ , |99|JT0T| that the proton diffuseness parameter in the WS model should be systematically 
increased in nuclei with large values of 2, in order to compensate for the reduction of the 
surface thickness due to the Coulomb potential. In particular, it has been shown in Refs. 



72| , |99H that the relative size of the Z—11A and 126 gaps is sensitive to the diffuseness in the 



proton average potential (the increased diffuseness gives rise to increased 2=126 shell gap). 
The difference between the HF and WS proton spectra seen in Figs. |] and |6| can thus be 
attributed to a non-self-consistent treatment of surface properties in the WS model. 

In the HFB+Gogny calculations of Ref. [9(| the Z—11A gap is small, around 1.4 MeV, 
and the 2=126 gap is even smaller (the spherical 2=138 gap is predicted to be around 
3 MeV). For the neutrons, the iV=184 and 228 shells are similar as in our calculations. 

Summarizing this section, our self-consistent calculations, based on the the Skyrme in- 
teractions SkP and SLy7, favor the 2=126 spherical shell closure over the 2=114 gap. This 
observation is consistent with the recent results of the experimental systematics of B(E2) 
rates ||102|| according to which 2=126 is presumably the dominant proton shell closure when 
the global properties of nuclei in the trans-actinide region are parametrized in terms of 
numbers of valence particles, which depend on the shell closure beyond the experimentally 
accessible range. 



IV. EQUILIBRIUM SHAPES 

It has been early noted |57| , |H2]JiTl3] , |104|| that the very existence of SHE elements is en- 
tirely due to shell effects. Indeed, for 2>100, the macroscopic energy becomes unstable to 
shape deformations due to the Coulomb repulsion. This instability is illustrated in Fig. [7] 
which displays -Emacro calculated using the diffused-surface Yukawa-plus-exponential mass 
formula (Sec. [II B|) as a function of /?2 for three elements with Z=108, 114, and 126. It 



is interesting to note that while for 2=108 and 114 the macroscopic energy is unstable 
to prolate distortions, for still heavier elements there develops an instability with respect 
to oblate deformations. (In reality, this apparent minimum is a saddle point, unstable 
to triaxial distortions [|105| , |106|| .) This softness against oblate deformations explains large 



oblate ground-state deformations in nuclei with 2>120 obtained in some calculations (see 



below). As shown in Refs. [11,106], similar results are obtained with the standard sharp- 
surface liquid-drop formula [ 107 1 ; the latter one is more rigid to deformation and gives rise 
to slightly wider and higher fission barriers. 

Since the macroscopic energy of the SHE is unstable with respect to spontaneous fission, 
it is the magnitude of the shell correction that governs their lifetimes. Moreover, since 
the macroscopic energy does not favor spherical shapes (contrary to the situation in lighter 
elements), the deformed shell effects play a very important role in the description of SHE. 
Practically any local deformed gap around the Fermi level can give rise to a local minimum 
in the potential energy surface (PES), i.e., to a multitude of coexisting minima. The detailed 
energy balance between the local minima (i.e., whether a nucleus is spherical, prolate, or 
oblate) is dictated by the distribution of spherical single-particle orbitals (see Sec. JT|). To 
illustrate this point, potential energy surfaces characteristic of nuclei from this mass region 
are shown in Figs. § (WS, 310 126i 84 ) and g (HF+SLy7 , 288 112 176 ). As one can see, in both 
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cases the pattern of PES is very irregular; it reflects shell fluctuations in E m i cro . According to 
HF+SLy7, the ground state minimum of 288 112 176 corresponds to a deformed prolate shape 
with <52o=15b (/32~0.11). There also appear a shallow oblate minimum with <52o~ — 25 b 
which, as we shall see below, corresponds to a saddle point (i.e., maximum in the 7-direction 
|62| , |106|| ), and a low-lying superdeformed (SD) minimum at Q2o~80b. The nucleus 310 126i84 
is predicted to be oblate (fa—— 0.2) by the WS calculation; the spherical minimum is about 
1 MeV higher in energy. Two excited local oblate pockets appear at fa=—0A and —0.6 in the 
region of very flat behavior of -E maC ro on the oblate side. There is also a local well-deformed 
prolate pocket (fa=0A5) beyond the first barrier on the way to fission. A similar complex 



pattern of PESs for SHE has been obtained in previous WS calculations [11] and in the 



recent HFB+Gogny calculations of Ref. [96 



A. Axial shapes 

In self-consistent calculations, the shape of a nucleus is characterized by expectation 
values of multipole operators rather than by shape deformations. In our axial HF calculations 
the expectation values of multipole moments 

Q X0 = 2r x P x (cose) (6) 

are determined for each self-consistent solution. For the sake of comparing the HF and WS 
calculations, one can introduce the equivalent deformation parameters, i.e., such that the 
uniform distribution of matter inside the nuclear surface defined by Eq. (|j) reproduces the 
HF expectation values Q\o=(Qxo) of multipole moments (||). By keeping the second-order 
terms in fa and fa, the expressions for the mass quadrupole and hexadecapole moments, 
which define the equivalent deformation parameters, are (cf. Ref. ||108|| ) 

where Rq = 1.2A 1 / 3 . In a similar way one can define the relation between neutron (proton) 
moments and the corresponding deformation parameters. 

Figure pi] displays, in the form of contour maps, the ground-state equilibrium values of 
the equivalent fa obtained in the HF+SkP and HF+SLy7 variants. The general pattern is 
very similar for both Skyrme interactions. Namely, for A<176 all the nuclei considered are 
predicted to be prolate, and they are spherical for A>180. These two regions are separated 
by a narrow band of transitional nuclei with small equilibrium deformations. The largest 
deformations, around /3 2 =0.26, are obtained for the lightest nuclei considered, i.e., for Z<114 
and A<162. The stabilizing role of deformation in the region of the known heaviest elements 
is well known p|-4T^[; it has been attributed to the deformed shell gaps at A=162 and Z=108. 



Similar ground-state deformations of the heaviest elements with Z<114 have recently been 



obtained in the HF+SIII calculations 109 
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Figure [11] shows the ground-state equilibrium values of the equivalent (3± deformations. 
Again, the isotopic and isotonic dependencies of f3± are very similar for both effective in- 
teractions. All nuclei which are deformed in the ground state have negative values of /J*. 
The largest hexadecapole distortions, /3 4 ~— 0.10, are calculated in nuclei around Z=112 and 
iV=168. 

The ground-state quadrupole deformations obtained in the WS model are displayed in 
Fig. O. For Z<118, the results are close to those calculated in the self-consistent models. 



The significant differences are obtained for nuclei around Z=122 and iV=166 (WS predicts 
very deformed oblate shapes with /32~— 0.4), and for the isotopes with Z>12A and those 
around Z=122 and iV=190 where WS yields deformed oblate shapes with /5 2 ~ — 0.2. In most 
cases, our present WS results are very close to the WS results of Ref. fll| (cf. Figs. 2 and 
3 therein). 

In order to relate our results to previous work, we present in Table |V| the quadrupole 
and hexadecapole ground-state deformations obtained in the WS model, finite-range 
droplet model (FRDM) of Ref. p7|| , extended Thomas- Fermi with Strutinsky-integral model 
(ETFSI) [ |110|| , and in our HF+SLy7 and HF+SkP calculations. For the lightest systems, 
264 108 and 270 110, all calculations agree with one another. The results of both HF variants 
are very close for all nuclei considered. 

As illustrated in Figs. § and [| and discussed above, one expects competition between 
near-lying prolate, oblate, and SD shapes. Hence, depending on those model details which 
can influence the magnitude of shell correction, the calculated ground state of a superheavy 
nucleus can be prolate, oblate, or SD. This is clearly seen in Table |V|, especially for the 
heavier isotopes with Z>110. 

In the FRDM calculations, the oblate minimum is favored for weakly deformed nuclei 
with Nml74. A similar result was obtained in a recent HFB+Gogny study of Ref. EH]. A 



SD ground state (/?2^0.5) is predicted in the FRDM for the very heavy systems with Z>118 
and iV>194. In the ETFSI model, many nuclei (practically all isotopes with Z>112) have 
very elongated ground states with /3 2 ~0.4 (cf. Table |V|). The HFB calculations of Ref. |96| 
are in good agreement with our results. Except for the few isotopes that are predicted to be 
oblate in the HFB+Gogny calculations, the general features of the potential energy surfaces 
are similar. 

The transition from deformed to spherical shapes predicted for iV>180 can be associated 
with the spherical neutron shell closure at iV=184 (see Sec. |T|). The effect of the Z=126 
spherical proton gap appears to be weakened as compared to the known nuclei. The reason 
for this is twofold. Firstly, as discussed in Sec. |HI| above, the gap itself is significantly 
quenched in the SHE. Secondly, unlike in the lighter elements, spherical shell effects are not 
further stabilized by the liquid-drop energy contribution; i.e., for Z=126 the macroscopic 
energy has a maximum at the spherical shape. Consequently, the spherical proton gap at 
Z=126 cannot be considered as truly magic: the neutron-deficient Z=126 isotopes are all 
predicted to be deformed. 



B. Triaxial shapes 

In order to illuminate the role of triaxial degrees of freedom, in selected cases we per- 
formed the constrained HF calculations in the (Q20) Q22) plane. Figures ^ and 14 display 
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the self-consistent PES's calculated for the HF+SLy7 interaction and for 288 112 and 310 126, 
respectively. Here, the triaxial deformation 7 is defined as: 

tan 7 = V2^. (8) 



As discussed in Sec. [IV A|, the nucleus 288 112 is weakly prolate-deformed in its ground 



state. The secondary oblate minimum seen in Fig. || corresponds to a saddle point. The 
first fission barrier represents an axial saddle point at Q 2 o=4:2b. A secondary SD minimum, 
at Q2o=80b, can be seen at the edge of the map and is more clearly displayed in Fig. |]. 
In short, the static properties of the nucleus 288 112 can be well described by calculations 
constrained to axial shapes. 

The situation is rather different in 184 126. For both interactions, the ground state is 
spherical and there is no indication for the oblate minima predicted in WS. The static path 
to fission goes through triaxial saddle points. There is only one saddle point at 7^— 30° for 
SkP and two triaxial saddle points and a secondary triaxial minimum for SLy7. Clearly, 
the microscopic study of the fission process in 184 126 would require a dynamical calculation, 
fully exploring the triaxial degrees of freedom. 



V. SHELL CORRECTIONS 

To compare relative strengths of shell effects in different models, it is convenient to 
introduce the shell energy, i? S heii- 111 WS, the shell energy is defined as the difference between 
the total binding energy and the spherical part of the macroscopic energy 

Esheii = E — £ macro (spherical) = E micTO (spherical) + E de{ , (9) 

i.e., it is the sum of the spherical microscopic energy and the total deformation energy E^ef- 
In order to extract -Egheii from the HF calculations, we applied a very similar procedure. That 
is, from the self-consistent ground-state energy we substracted the same -Emacro (spherical) 
as in the WS calculations. In addition, the obtained HF values of -E s heii were shifted by 
a (small) constant in such a way that their values at the Z=108, iV=162 nucleus become 
equal to the WS value. This procedure aims at comparing the shell effects calculated in 
different models irrespectively of a generic smooth Z and iV dependence and a constant 
offset. (Another way of extracting ^heii from the HF calculations would be to perform 
the Strutinsky renormalization procedure based on the HF single-particle energies ||111| . 



However, because of the unphysical description of the positive-energy discretized continuum, 
especially for the proton-rich SHE, the standard renormalization technique is not very useful 
|112|| and we did not pursue this approach further.) 



The contour maps of -Egheii are displayed in Fig. [T5|. Both HF calculations clearly show 
the presence of strong shell stability around the "doubly-magic" nucleus 310 126. This result 
markedly differs from both that of the WS variant, where the island of shell stability is 
concentrated around the nucleus 296 114, and the result of the FY model calculations by 



Moller and Nix p6 who obtained the maximum shell energy at 114. 



The position of the island of shell stability does influence the mass differences that 



largely determine proton and neutron separation energies (Sec. |VTD and particle emission 
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probabilities, and hence the partial lifetimes (see Sec. |VII| ). The fact that the values of 
£7 s heii have very different ranges in HF and WS calculations (Fig. 15) implies that the Z 



and N dependence of the total energy also differs; although both WS and HF calculations 
reproduce properties of known heavy elements, they extrapolate very differently with the 
particle numbers. 

VI. SEPARATION AND a-DECAY ENERGIES 

Based on the calculated values of binding energies (or mass excesses, M exc ), separation 
and nuclear reaction energies can be extracted. These quantities involve masses of two or 
more nuclei and they are sensitive probes of symmetry properties of the effective interaction. 

A. Two-neutron and two-proton separation energies 

Firstly, we calculate two-neutron and two-proton separation energies of even-even SHE. 
These quantities, defined in terms of M exc , 

S 2n = ~M exc (Z, N) + M CXC (Z, N-2) + 2M CXC (0, 1), (10) 
S 2p = ~M exc (Z, N) + M CXC (Z -2,N) + 2M exc (l, 0) (11) 

[M exc (l,0)=8.071MeV, M exc (0, 1)=7.289 MeV], probe the variations in the mass surface in 
the directions of N and Z. 

The calculated (HF and WS) two-proton and two-neutron separation energies are plotted 



in Figs. [16] and [L7], respectively. In spite of rather different values of -E s hcii obtained in the 
self-consistent calculations and in the macroscopic-microscopic approach, the behavior of 
two-particle separation energies is strikingly similar. The excellent agreement obtained for 
the known heaviest elements seems to hold for very exotic superheavy systems as well. 

The two-proton drip line, S2 P =0, (indicated by a thick solid line in Fig. [16]) goes very 
smoothly between the nuclei 264 110i54 and 314 126i8g. 111 agreement with the results presented 
in Sec. [TTT], this regular behavior of S 2p as a function of N does not indicate any strong shell 
effects. Interestingly, the nucleus 310 126i84, predicted to be strongly stabilized by the shell 
effects in HF, is calculated to be proton- unstable (see Fig. |16|). However, due to very large 
Coulomb barriers, nuclei predicted to have proton separation energies S p > —1.5 MeV are 
practically stable to proton emission, i.e., this decay channel can be practically disregarded 
as compared to other decay modes (cf. discussion in Sec. [VI I| ). A similar conclusion was 
drawn by Tondeur [92] who predicted proton instability and a long proton emission lifetime 
of 310 126 but estimated very short a- and /3 + -decay half-lives. 

As seen in Fig. [T7|, the isotopes considered are very far from the neutron drip line. A 
strong spherical shell effect at iV=184 is clearly visible, especially in HF calculations. The 
effect of the deformed closure at iV=162, particularly pronounced in SLy7, is also seen. 

B. Alpha-decay energies 

As mentioned earlier, the main decay mode for the known heaviest elements is a-decay. 
The calculated Q^-values, 
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Qa = M exc (Z, N) - M exc (Z - 2, N - 2) - M exc (2, 2) (12) 



[M exc (2, 2)=2.425 MeV], are displayed in Fig. |18[ The landscapes of Q a obtained in the 
HF+SkP and HF+SLy7 models are very similar. For the nuclei around Z—W2 and iV=162, 
the WS variant yields Q a -values extremely close to those of HF+SkP. This is illustrated in 
Table |VJ which displays calculated values of Q a for the known a-decay chain of 264 Hs and 
compares them to the data. The obtained good agreement with experiment and with the 
results of the WS model (which was optimized locally in this region) is a necessary condition 
for trusting predictions of our self-consistent calculations in the region of unknown SHE. 

A noticeable difference is obtained in the region around (Z=116, iV=184) where WS 
predicts values which are ~2 MeV larger than those in self-consistent calculations. This shift 
in Q a results from different predictions regarding proton shell effect (Z=114 vs. Z=126, 
see Sec. [TTT]) and has a direct consequence for a-decay half-lives in this region (see Sec. |VT l 
below) . 



C. Beta-decay energies 



In the region of /^-stability, both Q@- and Qec are negative. Since our calculations are 
limited to the even-even nuclei, in order to estimate the beta-decay energies 

Qp- = M exc (Z, N) - M CXC (Z + 1,N-1), 
Qec = M exc (Z, N) - M exc (Z - 1, N + 1), 



Q EC - 2m e 



(13) 



masses of the odd-odd systems have been obtained in the quasi-particle approximation 
|113| , |114|| . The final result for the beta-decay energy of the parent nucleus (Z, N) can be 
expressed in terms of binding energies of neighboring even-even nuclei and the "effective" 
LN pairing gaps: 



Qt 



AM H _ n - Ap (Z, N — 2) — X 2tP (Z, N — 2) — A^ N (Z, N - 2) - \ 2>n {Z, N -2) 



+ -[3B(Z,N)-B(Z,N 



Q 



EC 



-AM, 



2) - B(Z + 2, N) - B(Z + 2, N - 2)] , 
„_„ - A"f{Z -2,N)- \ 2 , P {Z -2,N)- A]f(Z -2,N)- A 2 , n (Z - 2, N) 

B(Z,N + 2)}, (14) 



+ - [3B(Z, N) - B(Z -2,N)- B(Z — 2, iV + 2 

where AM H _ n = m n — M H , M H and m n denote the hydrogen atom and neutron masses, 
respectively, and B is the binding energy of an even-even nucleus. 



The predicted /3-stability valley is shown in Fig. [19]. According to our calculations, the 
spherical nuclei from the neighborhood of 298 114 are beta-stable. On the other hand, the 
predicted "doubly- magic" nucleus 310 126 lies outside the stability valley. This result agrees 
with the conclusions of previous studies [^4| , p2| , p6| where the nuclei around 310 126 have also 
been calculated to be f3 + /EC-unstable. According to the estimates of Ref. fl4| , Gamow- 
Teller transitions are of the order of 0.3 s for Z w 128 and iV w 188. Kumar PDA applied the 



semi-empirical estimate for /3-decay and EC half-lives ||115|| and obtained the value of 18 s for 
310 126. Thus, in this region, the /? + /EC process is expected to be much slower as compared 
to the a-decay |)2] (see Sec. |VI1|) . On the other hand, for the neutron-rich isotones with 
iV^184 the j3~ emission is expected to compete with a-decay and fission [92]. 
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VII. PARTICLE DECAY LIFETIMES 



A. Alpha decay half-lives 

In order to estimate the a decay half-lives, T a , we have used the phenomenological 
formula of Viola and Seaborg [ |116| ] with the parameters of Ref . [ 25 1 , 



log 10 Tj s (sec) = (1.66175 Z p - 8.5166)/^^ - (0.20228 Z p + 33.9069), (15) 

where Z p is the atomic number of a parent nucleus and Q a is in MeV. The a-decay half-lives 
calculated in the HF+SkP, HF+SLy7, and WS models are shown in Fig. |20]for the even-even 
SHE. 

For the known heaviest elements, there is fairly good agreement between all three variants 
of calculation. In particular, our predictions for nuclei ZpsIIO-114 and iVfsl54-170 agree 
with those of Refs. [ |117| | and [Q. Especially pronounced in the HF+SLy7 model, there 
appears a region of increased stability associated with the deformed shell closure at iV=162 
|5|JTT[| . Table [VII] displays calculated values of T a for the nuclei which are the most promising 



candidates to be synthesized in the near future. 

The HF+SLy7 model predicts an increased stability with respect to the a-decay also for 
nuclei with iV"psl72 and Z ranging from 110 to 118. It originates from the deformed neutron 
gap that opens up for iV=172-174 at large negative values of (5^. The resulting rapid change 
of the slope of the shell energy as a function of particle number around N ^170 is seen in 
Fig- 

Both HF calculations differ significantly from the WS calculations for the heavier iV>180 
isotones and predict a considerably smoother decrease of T a with increasing proton number. 
This discrepancy can be attributed to the difference in the location of the spherical proton 
shell which is moved up from Z=114 to Z=12Q. The maximum in T a at iV^182 can be 
attributed to the spherical neutron shell effect at iV=184. Because the proton shell effect 
appears at higher values of Z, HF calculations predict longer alpha-decay half-lives for all 
nuclei with iV^182 as compared to the present WS calculations and the previous calculations 
of Refs. [ |iT| , P^[2"6T , |6"8"||l 1 7| . |1 1 5j based on the macroscopic-microscopic approach. On the other 
hand we note a similarity between our HF+Sly7 results and those of Ref. |96| within the 
HFB-Gogny model. 

It should be emphasized that Eq. (|T5|), based on the WKB approximation, provides a 



rather crude estimate of T a since it disregards many structure effects [|119|| such as deforma- 
tion, configuration changes, and so on. However, based on the experience with the known 
heaviest elements, the phenomenological expression (|15|) usually gives a right order of mag- 
nitude for T a ; for more precise estimates the microscopic calculations of a decay are called 
for. 



B. Proton emission half-lives 

Nuclei with negative proton separation energies are unstable to proton emission. The 
proton emission half-lives depend primarily on the proton separation energy and orbital 
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angular momentum but very weakly on the details of the parameters of the corresponding 
proton optical model potential [|120H122 . 



The proton decay lifetimes have been estimated using the WKB expression for the partial 
width: 



h 2 

r = Q 2 M— exp 
4// 



/•Tout 

2 / drk(r) 

Jr\ n 



(16) 



where 9 2 is the proton spectroscopic factor, // is the reduced mass, r in and r out are the 
classical inner and outer turning points, respectively, Af is the normalization factor, and 
k(r) is the wave number given by 



Kr) = ^\Q P -V p (r)\. (17) 

In Eq. (PI), Q p =— S p >0, and V p (r) is the average proton potential (including the Coulomb 
potential). 

In our calculations, the proton potential was approximated by that of the WS model, 
and it has been assumed that 8 2 =1 and the proton orbital angular momentum, £, is zero. 
Because of the neglect of the centrifugal barrier, the calculated proton half-lives, T p , are, 
in general, grossly underestimated. That is, the predicted values of T p should be viewed as 
lower limits of proton half-lives. 

The calculations have been performed for the even-even and odd-even nuclei with S p <0 
(see Fig. [TBp. Since the proton separation energy involves the binding energy of the odd- 
even system (Eq. ( PH) ) to calculate the latter quantity, we employed the quasi-particle 
approximation (see discussion before Eq. ([14]) ). The final result for the proton separation 
energy of the parent nucleus (Z, N) can be expressed in terms of the proton Fermi energy, 
X p , and the "effective" LN pairing gaps: 

= <\p(Z,N)-A p LN (Z,N)-\ 2p (Z,N) if Z is even 
^ p \x p (Z,N) + A p LN (Z,N) + X 2 p(Z,N) if Z is odd 1 ' 

Here the tilde denotes the average value, e.g., 

X P (Z, N) = i [X P (Z -2,N) + X P (Z, N)) . (19) 

The obtained proton-emission half-lives are displayed in Fig. [H]. For a given isotonic 
chain, the values of T p for odd-Z nuclei are dramatically lowered with respect to the even-Z 
neighbors. As seen in Eq. ([Tj|), this odd-even effect is due to the proton pairing energy. 

Due to large Coulomb barriers in SHE, nuclei with Q p <1.5MeV are practically proton- 



stable. As seen in Figs. 20 and 21, the a-decay always wins for all nuclei considered. That 



is, the contribution of proton emission to the total decay width can be practically neglected. 



VIII. FISSION BARRIERS 

The self-consistent static paths to fission have been computed in a few selected cases. 
The results for nuclei 288 112 and 310 126 are displayed in Figs. || and |22|, respectively. 
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The existence of a rather high fission barrier in 310 126 (about 14 MeV for SLy7 and 16 
MeV for SkP) is a consequence of increased stability of the Z=126 and iV=184 system. 
Interestingly, the patterns of PES obtained in both HF calculations for 310 126 are very 
similar. The increase in the barrier height for SkP is due to the lower energy of the ground- 



state configuration in this model (see Fig. [15]). Calculations based on the macroscopic- 



microscopic method predict the fission barriers to be significantly smaller (cf., e.g., Ref. [|TT|j ), 
due to the absence of the spherical Z=126 shell closure in those models. 

Fission barriers are reduced by a few MeV if the non-axial degrees of freedom are taken 
into account (see Figs. 13 and |2"^) . Nevertheless, the resulting barriers are high enough to 
suppress the fission process as compared to the alpha emission. 

The heights of axial fission barriers computed in the HFB+Gogny model are very 
similar to our axial results (non-axiality has been ignored in their calculations) . They found 
fission half-lives to be of the order of 10 6 years for the iV=184 isotones (with barriers of the 
order of 14 MeV). A substantial decrease of T s f due to the non-axiality can be expected. (A 
rough estimate gives a decrease by almost 8 orders of magnitude for SLy7.) In this context, 
however, it should be emphasized that the height of the fission barrier obtained along the 
dynamic path is expected to be higher than that obtained along the static path, and the 
effect of non-axiality is reduced in dynamic calculations |T5| , |123| |. Rather high fission barriers 
in SHEs (about 14 MeV) are also predicted by Tondeur in Ref. [p2 |. 

The SLy7 static fission path for the nucleus 288 112 displayed in Fig. || shows the presence 
of the SD minimum which gives rise to a two-humped structure. (Note that the non-axial 
degrees of freedom do not influence the height of barriers in this region (see Fig. [13]).) Both 
barriers (inner and outer) have approximately the same height of about 6 MeV. Because 
of the strong Coulomb repulsion between the fission fragments, the barrier of 6 MeV is the 
minimum to ensure the stability of a SHE. Hence the fission process will probably dominate 



in this region. A similar pattern for the static fission path has been obtained in Refs. [15,96 



In Ref. [IT], the inner barrier has been found to be lower by 2 MeV compared to our result. 



Such a decrease can be attributed to the weaker shell effects predicted in |LjJ . 

The self-consistent barrier calculations performed here for a few nuclei confirm the general 



pattern predicted in Refs. |IIJT5|. Namely, in the region of the deformed iV=162 neutron 



shell, the a-decay mode is favored due to strong shell effects which hinder the fission process. 
On the other hand, for nuclei with iV>166 and Z^\Y1 the fission mode dominates. This can 
be justified, e.g., by a very steep increase of the shell energy in the HF+SLy7 calculations 
leading to a rapid decrease of the fission barriers. Since our self-consistent models do not 
predict an increased stability in the neighborhood of Z=114 and A r =184, there will probably 
be strong competition between spontaneous fission and a-decay. Heavier isotones of iV=184 
are expected to decay mainly via the a emission as one is approaching the spherical magic 
shell Z=126. 

The SHEs with Z>110 were also analyzed by Kumar []30] within the dynamic deformation 
model. He obtained a-decay half-lives much larger than ours. For instance, his values of 
T a for 298 116 and 310 126 are 10 L85 y and 10 10 ' 04 y, respectively. According to this model, the 
nucleus 298 114 decays by fission and 298 116 decays by /3 + /EC. 
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IX. CONCLUSIONS 



We have applied the HF method with two different Skyrme interactions, SkP and SLy7, 
and the WS macroscopic-microscopic approach to the study of SHEs with 108<Z<128 and 
150<iV<192. The overall agreement with existing experimental data for the heaviest known 
elements is satisfactory. In particular, the HF calculations predict the increased stability 
due to the deformed shell effect at iV=162. However, the self-consistent and macroscopic- 
microscopic approaches give different results when extrapolated to still heavier elements. 
The fact that the former methods treat correctly the interplay between the nuclear and 
Coulomb properties, which is crucial for the physics of SHEs, gives us more confidence in 
self-consistent results. 

Predictions have been made for binding energies, shell energies, deformations, reaction 
energies, and half-lives. In general, both SkP and SLy7 Skyrme parametrizations give similar 
results. 

According to our study, the doubly-magic superheavy element is that with Z=126 and 
A=184. This system is unstable with respect to the /3 + /EC, proton, and a decays, and 
it fissions. However, its half-life is governed by the a-decay (T a <l //s). The alleged magic 
gap at Z—114 does not show up in our HF calculations. Its presence in the models, based 
on the macroscopic-microscopic approach, is due to the lack of a self-consistent treatment 
of surface properties (in particular: proton diff useness) . In the self-consistent theory we 
obtain an increased shell stability of spherical nuclei with N^184 and Z>114, which gives 
rise to increased a-decay half-lives and opens up a possibility to access experimentally new 
higher- Z superheavy elements. 
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TABLES 



TABLE I. Nuclear matter characteristics of the SkP and SLy7 Skyrme parametrizations. Stan- 
dard definitions of the energy per particle E/A, Fermi momentum kp, incompressibility K^, 
isoscalar effective mass m*/m, dipole enhancement factor k, and symmetry energy a sym are col- 
lected, e.g., in Ref. [35|. The surface energy a sur f is calculated using the expression given in Ref. 



E/A 


kp 




m* /m 


K 


Osym 


^surf 


(MeV) 


(fm" 1 ) 


(MeV) 






(MeV) 


(MeV) 


SkP -15.95 


1.34 


201 


1.00 


0.35 


30 


21.8 


SLy7 -15.91 


1.33 


230 


0.69 


0.25 


32 


21.5 



TABLE II. Strengths of the zero-range pairing interaction, Eq. (|T|), adjusted to reproduce the 
proton and neutron average pairing gaps in 254 Fm in the SkP and SLy7 variants of calculation. The 
resulting "effective" odd-even mass difference, A LN + A2, is compared to the average gap, Eq. (|2|), 
and the experimental value extracted from the odd-even mass differences [Q. All quantities are 



in MeV. 







SkP 




SLy7 








V 


A LN + A 2 


V 


A LN + A 2 


A 


^exp 


protons 


245 


0.823 


290 


0.817 


0.820 


0.768 


neutrons 


190 


0.645 


250 


0.642 


0.647 


0.645 



TABLE III. Odd-even mass differences A + A 2 calculated in the HF+SkP+LN approach 
compared to the average gaps, Eq. (||), and the experimental odd-even mass differences pi] for 
several actinide nuclei. All quantities are in MeV. 



neutrons protons 



Nucleus 


A LN + A 2 


A 


^exp 


A LN + A 2 


A 


^cxp 


232 Xh 


0.65 


0.68 


0.75 


0.89 


0.86 


1.03 


236-q 


0.65 


0.67 


0.64 


0.90 


0.85 


0.92 


240p u 


0.65 


0.67 


0.57 


0.85 


0.84 


0.76 


246 Cm 


0.70 


0.66 


0.58 


0.85 


0.83 


0.85 


250 Cf 


0.67 


0.65 


0.64 


0.84 


0.83 


0.78 
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TABLE IV. Size of the Z=IU and N=184 gaps (in MeV) obtained for the nucleus 298 114 
within various models. 



Model 



114 



184 



HF+SLy7 a 
HF+SkP a 
WS a 
WS [§6| 

HF+SIII |||l| 



HF+SKa Q 
BHF 
BHF M 



ETF [93 1 



RMF-NL1 |4| 
RMF H 



DDRMF H 



1.55 
0.63 
2.05 
2.5 
1.93 
2.2 
1.1 
1.13 
1.4 
1.47 
2.3 
1.8 



2.92 
2.36 
2.15 
2.14 
1.6 
2.7 
2.5 
2.64 
2.3 

3.0 
1.0 



a This work. 
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TABLE V. Ground-state deformations /?2 


(upper numb 


3r) and f3± (lower number" 


for selected 


superheavy nuclei predicted 


in the WS, FRDM, ETFSI, HF+SLy7, and HF+SkP models. 


Nucleus 


WS a 


FRDM b 


ETFSI 


HF+SLy7 a 


HF+SkP a 


264 108 156 


0.24 


0.23 


0.26 


0.26 


0.27 




-0.04 


-0.04 


-0.05 


-0.05 


-0.04 


270 110 160 


0.23 


0.22 


0.25 


0.25 


0.25 




-0.07 


-0.08 


-0.07 


-0.08 


-0.07 


278 H2 166 


0.20 


0.16 


0.43 


0.23 


0.22 




-0.10 


-0.06 


0.07 


-0.10 


-0.09 


288 H2 176 


0.09 


-0.06 


0.43 


0.11 


0.12 




-0.05 





0.05 


-0.07 


-0.05 


298 114 184 








-0.01 
















0.04 


-0.01 





294 120 174 


0.08 


0.08 


0.43 


0.09 


0.10 




-0.05 


-0.05 


0.05 


-0.05 


-0.05 


294 122 172 


-0.13 


-0.16 


0.43 


0.07 


0.09 




0.01 


0.03 


0.05 


-0.03 


-0.04 


306 126 180 


-0.24 


0.42 




0.01 


0.01 




0.04 


0.02 




-0.01 





310 126 184 


-0.21 















0.02 







-0.01 






a This work. 
b Ref. |7|. 
c Ref. BTTU]. 



TABLE VI. Calculated and experimental values of Q a for the 264 Mt— 


> 248 Fm a-decay chain. 


Nucleus 


HF+SkP 


HF+SLy7 


WS 


exp a 


264 Mt 


10.82 


10.81 


11.11 


10.80 


260 Nh 


10.16 


10.21 


10.12 


9.92 


256 Rf 


9.63 


8.77 


9.18 


8.95 


252 Nq 


9.04 


7.97 


8.51 


8.55 



*Ref. HI. 
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TABLE VII. Theoretical estimates for the T a values (in fis) for some nuclei 

Nucleus HF+SkP HF+SLy7 WS 

276 112 43 05 47 

280 114 1.3 3.8 1.4 

284 116 2.4 49.0 0.5 

288 118 2.8 8.6 0.9 
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FIGURES 



FIG. 1. Spherical single proton levels in Pb predicted in the WS model as a function of the 
Coulomb potential strength xqouI (for a definition see text). The single-particle levels are shifted 
by the energy Ae given by Eq. © with Vc O ui(r=0)=23.1597 MeV. 

FIG. 2. Same as Fig. | but for single proton levels in 208 126i 84 . Here Vfcoui(r=0)=31.2890 MeV. 

FIG. 3. Spherical neutron single particle levels in 298 114ig4 predicted in the HF+SkP, 
HF+SLy7, and WS models. 

FIG. 4. Same as Fig. |3| except for single proton levels in 298 114. 

FIG. 5. Same as Fig. || except for single neutron levels in 310 126. 

FIG. 6. Same as Fig. |3| except for single proton levels in 310 126. 

FIG. 7. The macroscopic energy -E maC ro (normalized to zero at spherical shape) given by the 
Yukawa-plus-exponential mass formula as a function of fa for 264 108, 298 114, and 310 126. At 
every value of fa, E ma _ cro is minimized with respect to deformations fa and fa. In addition to the 
instability to prolate distortions seen in all three nuclei, there also develops an instability for 310 126 
with respect to oblate deformations. 

FIG. 8. Potential energy curve for 310 126i84 obtained in the WS model as a function of fa. 
In order to emphasize shell effects, the spherical part of the macroscopic energy, -©macro, has been 
substr acted. At each value of fa, the total energy £7 to t (solid line) has been minimized with respect 
to higher- multipole deformations fa and fa. The macroscopic energy is also shown (dashed line). 

FIG. 9. Potential energy curve for 288 112i76 obtained with the Skyrme interaction SLy7 as a 
function of the quadrupole moment Q20- The total energy is normalized to zero at the spherical 
shape. Since the calculation was not restricted to axial shapes, the triaxial deformation 7 varies 
as a function of Q20 

FIG. 10. Contour map of the ground-state mass quadrupole deformation fa obtained with 
the Skyrme interactions SkP (top) and SLy7 (bottom), plotted as a function of Z and N. The 
deformation values at contour lines are indicated. 

FIG. 11. Contour map of the ground-state mass hexadecapole deformation fa obtained with 
the Skyrme interactions SkP (top) and SLy7 (bottom), plotted as a function of Z and N. The 
deformation values at contour lines are indicated. 
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FIG. 12. Ground-state mass quadrupole deformation obtained with the WS model plotted 
as a function of Z and N. 



FIG. 13. Potential energy surface of 288 112 in the (Q20, Q22)-plane calculated with the 
HF+SLy7 model. The difference between contour lines is 0.5 MeV. The local minima and sad- 
dle points are indicated by dots and crossed dots, respectively. 

FIG. 14. Same as in Fig. 13 except for the nucleus 310 126 calculated with the HF+SkP (top) 
and HF+SLy7 (bottom) models. 

FIG. 15. Shell energy Eq. (|), for SHE obtained with HF+SkP (top), HF+SLy7 (middle), 
and WS (bottom) model. Note different energy scale in the HF and WS calculations. 

FIG. 16. Two proton separation energies for even-even SHE obtained with the HF+SkP (top), 
HF+SLy7 (middle) and WS (bottom) models. The difference between contour lines is 1 MeV. The 
dashed lines indicate the nuclei with Sip < unstable to the two-proton emission. The positions 
of one-proton and two-proton drip lines are indicated by thick solid and dotted lines, respectively. 

FIG. 17. Two neutron separation energies for even-even SHE obtained with the HF+SkP (top), 
HF+SLy7 (middle) and WS (bottom) models. The difference between contour lines is 1 MeV. 



FIG. 18. Q a values for even-even SHE isotopes obtained with the HF+SkP (top), HF+SLy7 
(middle) and WS (bottom) models. The difference between contour lines is 0.5 MeV. Solid lines 
indicate nuclei with Q a <12.5MeV, i.e., those with alpha lifetimes greater than 1 [is (see Sec. VH 
and Fig. M). 



FIG. 19. The /3-stability valley calculated with the SkP (top) and the SLy7 (bottom) forces. 

FIG. 20. Alpha decay half-lives, T Q , for even-even SHE isotopes obtained with the HF+SkP, 
HF+SLy7, and WS models. 

FIG. 21. One-proton emission half-lives for even-even (solid line) and odd-even (dashed 
line) SHE obtained in the WKB approximation using the Q p values from the HF+SkP (left) 
and HF+SLy7 (right) models and assuming £=0. 

FIG. 22. Potential energy curve for the nucleus 310 126 obtained in the HF+SkP (top) and 
HF+SLy7 (bottom) models as a function of the mass quadrupole moment Q20 with (dashed line) 
and without (solid line) assumption of the axial symmetry. The stronger shell effect at Z=126 
in SkP as compared to SLy7 lowers the gs configuration in the former model. Consequently, the 
fission barrier in SkP is higher. 
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